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ABSTRACT 

One proposed mechanism for heating the solar wind, from close to the sun to beyond ~ 10 AU, invokes low- 
frequency, oblique, Alfven-wave turbulence. Because small-scale oblique Alfven waves (kinetic Alfven waves) 
are compressive, the measured density fluctuations in the solar wind place an upper limit on the amplitude of 
kinetic Alfven waves and hence an upper limit on the rate at which the solar wind can be heated by low- 
frequency, Alfvenic turbulence. We evaluate this upper limit for both coronal holes at 5Rq and in the near- 
Earth solar wind. At both radii, the upper limit we find is consistent with models in which the solar wind is 
heated by low-frequency Alfvenic turbulence. At 1 AU, the upper limit on the turbulent heating rate derived 
from the measured density fluctuations is within a factor of 2 of the measured solar wind heating rate. Thus 
if low-frequency Alfvenic turbulence contributes to heating the near-Earth solar wind, kinetic Alfven waves 
must be one of the dominant sources of solar wind density fluctuations at frequencies ~ 1 Hz. We also present 
a simple argument for why density fluctuation measurements do appear to rule out models in which the solar 
wind is heated by non-turbulent high frequency waves "sweeping" through the ion-cyclotron resonance, but are 
compatible with heating by low-frequency Alfvenic turbulence. 
Subject headings: 



1. INTRODUCTION 

A number of observations indicate that the solar wind un- 
dergoes spatially extended heating. For example, in situ mea- 
surements from satellites such as Helios & Voyager show that 
electrons and protons have non-adiabatic temperature profiles 
at heliocentric distances r ~ 0.3 — 50 AU (e.g., Freeman 1988; 
Gazis et al. 1994; Richardson et al. 1995; Cranmer et al. 
2009). Similarly, remote UVCS observations detect large ion 
temperatures that increase with heliocentric distance within 
a few solar radii of the Sun (Kohl et al 1998, Antonucci 
et al 2000). The observed heating likely plays a critical role in 
accelerating the solar wind to supersonic and super- Alfvenic 
speeds and also significantly impacts the plasma properties in 
the near-Earth space environment. 4 

Several heating mechanisms have been proposed to account 
for these observations, including magnetic reconnection (Fisk 
2003; Fisk et al 2003), plasma instabilities driven by an elec- 
tron heat flux or cross-field currents (Markovskii & Hollweg 
2002a, 2002b; Markovskii 2004), non-turbulent waves at fre- 
quencies comparable to the proton cyclotron frequency Q; 
(Abraham-Shrauner & Feldman 1977, Hollweg & Turner 
1978, McKenzie et al 1979, Hollweg 1981, Tu & Marsch 
1997, Marsch & Tu 1997, Hollweg & Isenberg 2002), tur- 
bulent waves extending from low frequencies up ~ £2, (Isen- 
berg & Hollweg 1983, Tu et al 1984, Marsch 1991, Chan- 
dran 2005), and turbulence that fluctuates only on time scales 
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4 Some studies, however, dispute that the observed temperature profiles 
imply heating, claiming that many of the observed properties of the solar 
wind can be explained by strong heating localized at the coronal base and 
the collisionless propagation of non-thermal particle distributions from the 
corona out into the interplanetary medium (see, e.g., Scudder 1992ab). 



> Sir 1 . (Coleman 1968; Matthaeus et al 1999; Dmitruk 
et al 2001, 2002; Cranmer & van Ballegooijen 2005; Cran- 
mer, van Ballegooijen, & Edgar 2007; Verdini & Velli 2007; 
Chandran et al 2009). This paper focuses on this last class of 
models, those involving low-frequency turbulence. 

These models hypothesize that Alfven waves are launched 
into the corona by photospheric motions driven by solar con- 
vection (e.g., Cranmer & van Ballegooijen 2005). Some of 
the Alfven waves are reflected by the spatial gradient of the 
Alfven speed, and interactions between oppositely directed 
Alfven waves then cause the wave energy to cascade to small 
scales perpendicular to the magnetic field (X±). However, the 
correlation lengths of Alfven wave packets in the direction 
of the background magnetic field remain comparatively 
long, and the Alfven-wave frequencies remain <C £2,. When 
the wave energy cascades to perpendicular scales X± com- 
parable to the proton gyroradius p,, the wave energy begins 
to dissipate, heating the corona. The role of turbulence is to 
transform the initial wave energy into a form - fluctuations 
with very small perpendicular scales - that can be efficiently 
dissipated by kinetic processes. 

One prediction of these models is that the character of 
the turbulent fluctuations changes with decreasing X±. At 
X± 3> dj, the fluctuations are non-compressive Alfven waves, 
where di — Va/Q is the ion inertial length. 5 On the other 
hand, in low-p plasmas at p, < X± < di, the fluctuations are 
highly compressive, in the sense that 8«/«o > 8Z?/Z?o, where 
8n and no are the fluctuating and background electron densi- 
ties, and SB and Bo are the fluctuating and background mag- 
netic fields. For p ~ 1, 8«/«o ~ &B/#o at X± ~ p*. At even 
smaller scales, X± < p,-, the ion and electron velocities decou- 
ple for all p, the waves become dispersive, and the fluctuations 
produce parallel electric and magnetic field perturbations that 
cause the waves to damp. At X± < p,, these solutions to the 

5 For reference, dj ~ p _1 ' 2 p j , p = Sjtp/fi 2 <C 1 in the corona, and for 
typical coronal-hole parameters d, ■ ~ 3 x 10 5 cm at r ~ 5i? .-,; by contrast, 
p ~ 0.2 — 1 in the near-Earth solar wind at ~ 1 AU and d, ; ~ 3 x 10 6 cm. 
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Alfven-wave branch of the linear dispersion relation are called 
kinetic Alfven waves (KAWs). Many of the measured proper- 
ties of the magnetic and electric field fluctuations in the near- 
Earth solar wind are consistent with KAW turbulence on small 
scales (e.g., Howes et al. 2008a, 2008b; Sahraoui et al. 2009, 
and references therein). 

A very important, but somewhat indirect, test of the low- 
frequency turbulence model is provided by measurements of 
ion and electron temperatures, both in situ and in coronal 
holes (the open-field-line regions from which the fast so- 
lar wind is thought to emanate). Observations with UVCS 
show that 7j_ 3> 7]| for minor ions such as +5 in coronal 
holes, where 7j_ and 7j| are the temperatures corresponding 
to thermal motions perpendicular and parallel to the back- 
ground magnetic field Bo (Kohl et al 1998, Antonucci et al 
2000). It has not yet been possible to unequivocally deter- 
mine whether protons have the same temperature anisotropy 
in the corona. However, in situ measurements of the fast solar 
wind show that 7j_ > 7ji for the core of the proton distribution 
function, although the anisotropy is less pronounced for larger 
P (Marsch, Ao, & Tu 2004; Hellinger et al. 2006). 

It is not yet clear whether low-frequency Alfvenic turbu- 
lence can explain perpendicular ion heating and the prefer- 
ential heating of minor ions. Dissipation of Alfvenic turbu- 
lence occurs at small scales, at which the amplitude of the 
fluctuations is very small. The hot-plasma dispersion relation 
for small-amplitude waves thus provides a plausible first esti- 
mate for how KAW turbulence dissipates in the collisionless 
corona and solar wind. This predicts negligible damping for 
oblique Alfven waves with k±p, <C 1, where k± and kn are 
the wavevector components perpendicular and parallel to Bo 
(Barnes 1966). As a result, the energy cascades to scales < p,-. 
The damping of the resulting KAWs with k± ^> \k\\ | and fre- 
quencies <C £2, arises from Landau and/or transit-time damp- 
ing and leads to parallel electron heating for p < 1, not per- 
pendicular ion heating (Quataert 1998; Quataert & Gruzi- 
nov 1999; Leamon et al 1999; Cranmer & van Ballegooi- 
jen 2003; Gary & Nishimura 2004; Howes et al 2008a). Nev- 
ertheless, a number of mechanisms have been proposed that 
might produce perpendicular ion heating from low-frequency 
turbulence, including heating by reconnection electric fields 
(Dmitruk, Matthaeus, & Seenu 2004; see, however, Lehe, Par- 
rish, & Quataert 2009), secondary plasma instabilities trig- 
gered by the KAW velocity shear (Markovskii et al 2006), ion 
heating by electron phase-space holes generated by the heat- 
ing of electrons (Matthaeus et al 2003; Cranmer & van Balle- 
gooijen 2003), and perpendicular ion heating through stochas- 
tic ion orbits (Johnson & Cheng 2001). An outstanding prob- 
lem is whether any of these mechanisms (or others) can quan- 
titatively account for the measured perpendicular ion heating 
in the context of low-frequency turbulence models - this is 
a particularly important question given the growing body of 
evidence that the magnetic and electric field fluctuations at 
~ 1 AU are consistent with low-frequency anisotropic Alfven 
waves and KAWs (e.g., Sahraoui et al. 2009). 

The goal of this paper is to quantitatively assess a second 
test of the importance of low-frequency turbulence in the solar 
wind: the measured power spectrum of density fluctuations. 
Coles & Harmon (1989) measured the spectrum of electron 
density fluctuations in the corona at radii as small as 5Rq us- 
ing Venus as a background radio source. In addition, mea- 
surements of density fluctuations in the solar wind at ~ 1 AU 
have been carried out by a variety of methods (e.g., Celnikier 



et al. 1983, 1987; Hnat et al. 2005; Kellogg & Horbury 2005). 
Because oblique Alfven waves become increasingly compres- 
sive for k±di ~ 1, the observed density fluctuation spectrum 
constrains the spectrum of oblique Alfven waves and provides 
an upper limit on the heating rate e from Alfvenic turbulence. 
In § |2] we summarize the predicted density fluctuations in- 
duced by low-frequency Alfvenic turbulence, and describe 
some of the remaining uncertainties in these predictions (see 
also Schekochihin et al. 2009). We then show how the mea- 
sured density fluctuations in coronal holes and at ~ 1 AU con- 
strain the heating rate e due to Alfvenic turbulence ($3). In §|4] 
we describe the "sweeping" model for coronal heating (e.g., 
Schwartz et al 1981; Axford & McKenzie 1992; Marsch & 
Tu 1997; Ruzmaikin & Berger 1998), in which the corona is 
heated by cyclotron damping of kHz waves launched directly 
from the Sun; we present a simple explanation for why radio 
observations do not rule out turbulent heating models, even 
though they appear to rule out the sweeping model (as was 
shown by Hollweg 2000). In §[5]we summarize and discuss 
our results. Our analysis and results are broadly similar to 
those of Harmon & Coles (2005), but we consider a different 
model for the turbulent cascade of Alfven waves, and apply 
the density fluctuation constraint both to remote observations 
of turbulence in coronal holes and to in situ measurements of 
turbulence in the near-Earth solar wind. 

2. PREDICTED DENSITY FLUCTUATIONS IN LOW-FREQUENCY 
ALFVENIC TURBULENCE MODELS 

In this section we briefly summarize the density fluctuations 
produced by low-frequency Alfvenic turbulence (see Lith- 
wick & Goldreich 2001 and Schekochihin et al. 2009 for a 
more comprehensive discussion). We assume that the Alfven 
wave power spectrum follows the "critical balance" theory of 
Goldreich & Sridhar (1995). Initially we also assume that the 
turbulence is "balanced," i.e., that there are equal fluxes of 
Alfven waves propagating towards and away from the Sun in 
the plasma frame, or equivalently, that there is zero cross he- 
licity; we discuss the effects of a non-zero cross helicity at the 
end of the section. 

Although Alfven waves themselves are not compressive 
when k±dj ^> 1, low-frequency Alfvenic turbulence nonethe- 
less produces significant density fluctuations, via two differ- 
ent physical processes. First, as energy cascades to scales 
< dj, Alfven waves transition to KAWs, which are compres- 
sive (see Fig. [3] discussed below). Second, both slow waves 
and entropy modes are passively mixed by the Alfvenic cas- 
cade (Lithwick & Goldreich 2001). Because the Alfvenic 
fluctuations have an anisotropic Kolmogorov spectrum, the 
density fluctuations associated with the entropy modes and/or 
slow waves also have an anisotropic Kolmogorov spectrum. 

For the collisionless conditions appropriate to the solar 
corona and solar wind, the entropy modes and slow waves 
are both damped, but their damping rates are ~ knv t h,p, where 
v t h p is the proton thermal speed. 6 The cascade rate is ~ 
k\hvk ~ vpjc\\, where Sv^ is the rms amplitude of the velocity 
fluctuation at a perpendicular scale k~ l (where k ~ k^ since 
k± ^> feii). The cascade rate is thus comparable to the lin- 
ear Alfven wave frequency. For |3 < 1, the slow waves and 
entropy modes cascade faster than they are damped, and the 

6 Under these conditions, there are really two entropy modes, one associ- 
ated with the electrons and one associated with the protons. The former is 
strongly damped, at a rate ~ kv,i, e , while the latter is less strongly damped. 
We focus on the latter throughout. 
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the kinetic energy spectrum by specifying the ratio 



FIG. 1 . — One-dimensional density fluctuation power spectrum 5>i e pro- 
duced by balanced low-frequency Alfvenic turbulence in the near-Earth solar 
wind with p, = 0.43 and Tj/T e = 1. Total density fluctuation spectra (solid) 
are shown for / = 0. 1 , 1 , 10, with the separate passive-scalar component (dot- 
ted) and "active" KAW component (dashed) shown explicitly for / = 0.1. 



passive scalar contribution to the density fluctuations extends 
to small scales ~ p, (Lithwick & Goldreich 2001). We discuss 
the case of p > 1 below. 
Figure [TJ summarizes our prediction for the 1 -D density 

fluctuation spectrum Q>ne associated with balanced (zero 
cross-helicity) low-frequency Alfvenic turbulence. The rms 
density fluctuations associated with the slow waves and en- 
tropy modes are cascaded by the Alfvenic fluctuations like a 
passive scalar and, at perpendicular scale k , are denoted by 

8n[ ps '; the rms density fluctuations that actively arise from the 

kinetic Alfven wave compressions are denoted by 8n[ comp '. 
At large scales, the passive-scalar density fluctuation spec- 
trum is chosen to be proportional to the perpendicular kinetic 
energy spectrum 7 associated with the perpendicular velocity 
fluctuation 8v±/ ( . The latter is determined using the analytic 
cascade model of Howes et al. (2008a), which is based on the 
assumptions of local nonlinear energy transfer in wavenumber 
space together with the critical balance between linear propa- 
gation and nonlinear interaction times. The analytic model 
smoothly transitions from anisotropic Kolmogorov Alfven- 
wave turbulence at k±pi <C 1 to anisotropic KAW turbulence 
at k±pi ^> 1; it agrees well with numerical simulations of ki- 
netic turbulence in the (limited) comparisons available to date 
(Howes et al. 2008b). 

The relative magnitude of the passive and active density 
fluctuations is uncertain and may vary with position (and 
time) in the solar wind. We determine the constant of pro- 
portionality between the passive-scalar density spectrum and 

7 For small scales k^pt > 1, the passive-scalar density spectrum follows 
perpendicular magnetic energy spectrum rather than the perpendicular kinetic 
energy spectrum. 
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where 2n/ko is the driving scale or outer scale of the tur- 
bulence. Figure [TJ presents 1-D density spectra from solu- 
tions of the cascade model for near-Earth solar wind con- 
ditions. Plasma parameters for this figure have been cho- 
sen to correspond to period II of Celnikier et al. (1987): 
Bo = 1.5 x 10~ 4 G, m = 18 cm" 3 , T t = T e = 1.45 x 10 5 K, 
and v sw = 460 km/s, giving p,- = 2 x 10 6 cm. We assume a 
driving scale 2n/ko — 2x 10 cm, giving cascade model pa- 
rameters p, = 0.43, Ti/T e = 1, and /t p, = 6 x 10~ 5 . 

Figure [TJ shows the total density fluctuation spectra (solid) 
for / = 0.1,1,10, with the separate passive-scalar compo- 
nent (dotted) and "active" KAW component (dashed) shown 
explicitly for the / = 0.1 case. The / = 0.1 case demon- 
strates that the passively mixed density fluctuations have a 
Kolmogorov power-law spectrum at large scales (small k) 
with a break at k±pi ~ 1 where the turbulence transitions 
to dispersive KAWs. On small scales, however, the "active" 
density fluctuations due to KAWs become important and can 
dominate over the passive contribution. The density fluc- 
tuation spectrum associated with the KAWs is a factor of 
~ k\ flatter than Kolmogorov for k^pi < 1. Indeed, the ID 
density-fluctuation spectrum from KAW compressions is ris- 
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ing, oc k~{ " , tor &xP; ^ 1- This is because the density pertur- 
bation due to a linear KAW is °< k± times the velocity pertur- 
bation in this limit (see eq. J2) below). 

The / = 0.1 result in Figure [TJ may be compared directly 8 
to the measured density spectrum for this period of solar wind 
turbulence shown in Figure 5 of Celnikier et al. (1987). The 
predicted density spectrum, combining both a passive-scalar 
contribution and compressive KAWs, qualitatively reproduces 
the measured spectrum. Note that the assumption of balanced 
turbulence (zero cross helicity) is not an unreasonable as- 
sumption for this interval of relatively slow and dense solar 
wind. In addition, the value / = 0. 1 is consistent with several 
observations that are discussed following equation (01, below. 

In coronal holes, where p, <C 1, the compressive motions 
associated with the KAWs become comparatively stronger. 
Figure [2] presents cascade-model results analogous to Fig- 
ureQJfor parameters appropriate to coronal holes: P, = 0.01, 
Tj/T e = 2, and AioP; = 10~ 8 . 9 The absolute normalization of 
the density fluctuations was chosen to roughly match the ob- 
servations of Coles and Harmon (1989) in coronal holes. Total 
density fluctuation spectra (solid) are shown for / = 0. 1 , 1 , 10, 
with the decomposition into passive-scalar component (dot- 
ted) and "active" KAW component (dashed) included for each 
case. The results in Figure [2] are similar to those in Figure [TJ 
except that the active KAW component of the density fluc- 
tuations is more prominent at low P,-. The questionable as- 
sumption of balanced (zero cross-helicity) turbulence in coro- 
nal holes is discussed in detail at the end of this section. 

Figure [3] focuses on the density fluctuations due to linear 
KAWs, showing results as a function k±pj for several p,-. For 

8 Note that a value of k±pi = 1 in Figure [T]corresponds to a frequency of 
3.6 Hz in Figure 5 of Celnikier et al. (1987). 

9 In the cascade model, ko is not the actual outer-scale wavenumber. In- 
stead, in the case of the corona, it is the much smaller wavenumber at which 
the anisotropic power spectrum would extrapolate to isotropic fluctuations 
with SB = B (Howes et al 2008a). 
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FIG. 2. — One-dimensional density fluctuation power spectrum 5>i e pro- 
duced by low-frequency Alfvenic turbulence in coronal holes with P, = 0.01 
and Tj/Tg = 2. Total density fluctuation spectra (solid) are shown for / = 
0.1,1,10, with the separate passive-scalar component (dotted) and "active" 
KAW component (dashed) shown explicitly for the each case. 




FIG. 3. — Density fluctuations relative to the perpendicular electron velocity 
fluctuations (8n£/no)/(Svxfe/v>i) for (3; = 0.01,0.1,1,10 (solid) based on 
linear Vlasov-Maxwell theory for kinetic Alfven waves. Dotted lines are the 
predictions given by equation (2} with Y; = 1 . 



P < 1, the relation between the density and velocity fluctua- 
tions of a KAW at wavevector k, denoted 8«jt and Svjjt, can 
be derived analytically and is given by (Hollweg 1999, Eqn. 
51) 



8n k 



no 



k±di 



(1- 



-Jik 2 " 2 



Pi) 



5v 



Ik 



VA 



(2) 



assuming that k± ^> | and CO <C O,-, where CO is the wave fre- 
quency and Ji is the adiabatic index of the ions. In addition to 
showing the (normalized) ratio of density fluctuations to ve- 
locity fluctuations from linear kinetic theory (solid), Figure [3] 
also shows the values predicted by equation (O (dotted) with 
Yi = 1. This figure demonstrates that the density fluctuations 
are suppressed by a factor of ~ p when p > 1, as expected 
because the fluctuations become increasingly incompressible 



for high p. An additional change in the predicted density fluc- 
tuations for p > 1 - as can occur in the solar wind at ~ 1 AU - 
is that the damping time of slow waves and entropy modes be- 
comes shorter than their expected cascade times. Because the 
solar wind is collisionless, this occurs even on large scales. 
As a result, the passive contribution to the density fluctua- 
tions may be significantly diminished when p > 1 (Lithwick 
&Goldreich2001). 

Equation <j2j> can be used to estimate the ratio between 
the passive-scalar density fluctuation and the density fluc- 
tuation arising from KAW compressions. We assume that 



8« 



(comp) 



/Svjjt equals the ratio 8nk/8vxk for linear KAWs 



given by equation (|2). When P < 1 and the cross helicity is 
zero, Bn^ /Svjjt is independent of k within the inertial range. 
Assuming that Sn[ p /Svjjt remains roughly constant from the 
outer scale all the way to k± ~ pj~ , equation (01 implies that 
for balanced turbulence with p < 1 



(comp) 



5»[ ps) 



(P/) 



-1/2 
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k=P7 



At 1 AU, the value of Sn[ /«o at k = ko ranges from roughly 

0. 1 to 0.3, as the assumed outer-scale time scale (as measured 
in the spacecraft frame) of the turbulence is varied from 1 — 12 
hrs (Tu & Marsch 1995, Fig. 2-14). If we set v A = 77 km/s 
as in Celnikier's period II data and bv±/ ( ~ 30 km/s at k = ko 
at 1 AU (see Cranmer & van Ballegooijen 2005, Fig. 9), then 
/ ~ 0.07 - 0.6. Equation © thus suggests that KAW com- 
pressions usually dominate the density fluctuations on small 
scales (k ± pi ~ 1) at 1 AU. 

It is important to note several theoretical uncertainties in our 
predictions for the density fluctuations that should accompany 
low frequency Alfvenic turbulence in the solar wind. First, 
the relative amplitudes of the slow waves, entropy modes, and 
Alfven waves on large scales are not precisely known. There- 
fore, we present several values of / in Figures [TJ and [2] A 
second, and more significant, limitation is that we have as- 
sumed that the turbulence is "balanced" (zero cross-helicity), 

1. e., that the energy in waves propagating towards the sun (in 
the solar-wind frame) equals the energy in waves propagating 
away from the sun. This is generally not observed to be true 
at ~ 1 AU (particularly in the fast wind; Grappin et al. 1990; 
Tu & Marsch 1990). Closer to the Sun there is an even larger 
difference between the energies of Sunward and anti-Sunward 
waves (Roberts et al 1987; Bavassano et al 2000). Non-zero 
cross helicity or "imbalance" affects strong Alfvenic turbu- 
lence in several ways. For example, nonlinear interactions 
among Alfven waves occur only between waves propagating 
in opposite directions in the plasma frame (Iroshnikov 1963, 
Kraichnan 1965), and so if the energy in waves propagating 
towards the Sun is very small, then the energy cascade rate 
and turbulent heating rate also become small (Dobrowolny 
et al 1980, Hossain et al 1995, Chandran et al 2009). Fi- 
nite cross helicity may also modify the wavenumber scalings 
of the inertial-range power spectra of the density, magnetic 
field, and velocity. Because the "minority" Sunward Alfven 
waves and the passive scalar fluctuations are both cascaded 
by the "dominant" anti-Sunward Alfven waves, the passive 
scalar spectrum is expected to have the same inertial-range 
scaling as the Sunward waves (Lithwick & Goldreich 2003; 
Chandran 2008b). In some studies of Alfvenic turbulence 
with cross helicity (e.g., Grappin et al 1983; Chandran 2008a; 
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Beresnyak & Lazarian 2009) the minority Alfven waves have 
a shallower power spectrum than the dominant Alfven waves, 
suggesting that the passive scalar spectrum is shallower than 
the magnetic spectrum in highly imbalanced turbulence. This 
finding may be related to the shallow density spectra seen in 
radio observations of the corona, in which the turbulence is 
expected to be highly "imbalanced" (Cranmer & van Balle- 
gooijen 2005, Verdini & Velli 2007). 10 However, a number 
of other studies find that the inertial-range spectra of the mi- 
nority waves and dominant waves scale with wavenumber in 
the same way (Lithwick, Goldreich, & Sridhar 2007; Perez & 
Boldyrev 2009; Podesta & Bhattacharjee 2009). Moreover, 
the different studies cited above disagree over whether the 
spectra of the minority and dominant waves are equal at the 
dissipation scale ("pinning"). Because imbalanced Alfvenic 
turbulence is still not fully understood, and is likely the norm 
in the solar wind, we are only able to derive upper limits on 
the turbulent heating rates from the density observations, as 
discussed further below. In addition, this uncertainly implies 
that the precise power-law scalings for 4> n |? in Figures [T]and 
|2] should not be taken too literally, although we believe that 
our conclusions about the relative contribution of the active 
and passive density fluctuations are robust. 

3. OBSERVATIONAL CONSTRAINTS ON LOW-FREQUENCY 
ALFVENIC TURBULENCE MODELS 

In this section, we discuss observational constraints on den- 
sity fluctuations in the solar corona and solar wind, and their 
implications for low-frequency Alfvenic turbulence models. 

3.1. Coronal Holes 

Coles & Harmon (1989) analyzed the spectral broadening 
of Arecibo radar observations of Venus near superior conjunc- 
tion to determine the three-dimensional power spectrum of 
electron density fluctuations <P ne (r,k) (treated as an isotropic 
function of wave vector k) at a range of heliocentric dis- 
tances r in the slow solar wind. They were able to fit <b ne {r,k) 



at r = 5R & - 

An upper limit on the rms amplitude of the Alfvenic veloc- 
ity fluctuation at perpendicular scale fcr , denoted Svjj, can be 

obtained by assuming that the density fluctuations at scale k^ 1 
arise entirely from KAWs. Using the linear eigenfunctions of 
KAWs, we can write 



VA 



< 



no 



(7) 



If the compressibility of the KAWs noticeably affects the den- 
sity spectrum from 10" 7 cm" 1 < k < 10" 5 cm" 1 , as conjec- 
tured above and as is suggested by Figured then Sv^ may be 
close to the upper limit given in equation (0; this is because 
the fraction of <t>„ c that arises from KAWs increases at larger 
k (Fig. |2]). To evaluate the right-hand side of equation (0 we 
assume y, = 1 and 7} = 2.0 x 10 6 K and adopt the coronal-hole 
magnetic field model of Cranmer & van Ballegooijen (2005) 
(their eq. 2), which gives B = 5.8 x 10" 2 G at r — 5R Q . 
These parameters give p, = 2.3 x 10 4 cm, d\ = 3 X 10 5 cm, 
va = 1.6 x 10 8 cm/s, and 



tk<o.35 f5 ' lk ' 



"0 



or, equivalently, 



S Vjfe . <8.4km/s. 



(8) 



(9) 



At k — kj, the density spectrum at r — 5Rq is still close to the 
value obtained by extrapolating a power-law fit to <P ne for val- 
ues of k between 10" 7 cm" 1 and 10" 6 cm" 1 . We can thus as- 
sume that most of the cascade power is still present at k± — kj 
and that most of the dissipation occurs at k± > kj. More- 
over, because fc/p,- ~ 0.2, the kinetic energy and magnetic en- 
ergy of Alfvenic fluctuations at k± = kj are comparable, as 
in incompressible MHD, but not like the short-wavelength 
regime k±pi 3> 1, in which the magnetic energy dominates, 
at r = 5Rq with one power law at k < 10"' cm" 1 , a slightly The energy density of Alfvenic fluctuations at scale k~[ l is 



shallower power law at 10 7 cm 1 < k < kj, and an exponen- 
tial or Gaussian at k > ki, where 

h ~ 10" 5 cm" 1 (4) 
is the "inner-scale" wave number at r = 5Rr:,. They found that 
at k = k, and r = 5R Q , <i>„ e ~ 9.0 x 10 3 cm" 6 km 3 (see their 
Fig. 4). We focus on the inner scale for reasons that will 
become clearer below. Coles et al (1991) found that <$> ne is a 
factor of ~ 15 smaller in coronal holes than in the slow wind, 
and thus we set 

<S>ne(r = 5R e ,k= 10" 5 cm" 1 ) ~ 6.0 x 10 2 cm" 6 km 3 . (5) 
The rms electron density fluctuation Sn^. is given by Sn? ~ 

4nk?<E>„ e (ki), which implies 8n^ ~ 87 cm" 3 at r = 5R Q . We 
estimate the coronal-hole electron density from Eqn. (4) of 
Feldman et al (1997), which gives n e = 5.9 x 10 3 cm" 3 at 
r — 5Rq; this is very close to the value inferred by Fisher & 
Guhathakurta (1995) from observations taken with the Spar- 
tan 201-01 coronagraph. This value for the background den- 
sity then gives 

^1.5xl0" 2 (6) 
n 



thus ~ pSv^ . The time required for the fluctuation energy at 
k± = ki to cascade to k± > 2kj, denoted t c , satisfies the in- 
equality 



(10) 



The right-hand side of equation ( TTOb is the shearing time scale 
for an Alfvenic velocity fluctuation at k± = ki with rms ampli- 
tude 8v^.. This is a lower-limit on f c because in incompressible 
MHD, the anti-Sunward waves are sheared by waves propa- 
gating towards the Sun (Iroshnikov 1963; Kraichnan 1965), 
as summarized in ^2] Thus, if the Sunward waves are much 
less energetic than the anti-Sunward waves at k± = kj, the 
time required for Sunward waves to shear and substantially 
distort anti-Sunward waves is much greater than (fc/Sv^.)" 1 
(Lithwick, Goldreich, & Sridhar 2007; Beresnyak & Lazar- 
ian 2008; Chandran 2008a). The cascade power in low- 
frequency Alfven-wave turbulence at r — 5Rq, denoted e is 
roughly pv? /V c . Given equation ( flOl ), we can write 



e < copfe(8v|, 



(11) 



10 For example, Markovskii & Hollweg (2002) fit the one-dimensional 
density spectrum at r = 5R@ and at k < 10" 7 cm" 1 with a power law of 
the form &"' 2 and a power law of the form &" 2 ' 3 at 10" 7 cm"' < k < 
10" 5 cm" 1 . The density fluctuations at k < 10" 7 cm"' in their study are 
thought to con'espond to inertial-range passive-scalar fluctuations. 



where cq — 0.25 in strong incompressible MHD turbulence 
with equal fluxes of Alfven waves propagating parallel and 
anti-parallel to the background magnetic field (Howes et al. 
2008a). Substituting equation (O into equation (fTTT i. we find 



£< 1.5 x 10 



erg cm s 



(12) 
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This upper limit can be compared with the parameterized 
heating rates employed in empirical models of the fast so- 
lar wind. Allen et al (1998) constructed a series of two- 
fluid models with heating rates chosen to match in situ 
measurements of the fast wind at 1 AU as well as the 
coronal-hole density profile inferred from the brightness pro- 
file of electron-scattered, polarized white light (Fisher & 
Guhathakurta 1995). In these models, denoted SW2, SW3, 
and SW4, the total (electron plus proton) heating rates at r — 
5Rq were 3.1 x 10~ 9 ergcm~ 3 s , 1.4 x 10~ 8 ergcm~ 3 s _1 , 
and 6.8 x 10~ 9 ergcm~ 3 s _1 , respectively. Esser et al (1997) 
constructed similar models with higher temperatures at the 
coronal base and lower heating rates. The value of e at 
r = 5Rq was 2 x 10~ 10 ergcm~ 3 s~' in their model A and 
8 x 10~ 10 ergcm~ 3 s _1 in their model B. 11 Since these mod- 
els (by construction) provide a reasonable fit to the observed 
properties of the fast wind, we take the actual heating rate 
in coronal holes at r = 5Rq to be in the range spanned by 
these models, i.e., between 2 x 10~ 10 erg cm s -1 and 1.4 x 
10~ 8 ergcm _3 s _1 . Since the upper limit in equation ( flZb is 
above this range, we conclude that a model in which the fast 
wind is accelerated by heating from low-frequency Alfvenic 
turbulence is consistent with the radio observations. Future 
investigations could provide tighter constraints on the frac- 
tion of Sn*;. that arises directly from KAWs as well as the ratio 
of Sunward to anti-Sunward wave energy at k± = k,, which 
could in principle lower the upper limit on e that is implied by 
the density fluctuation measurements. 

3.2. The Solar Wind at ~ 1 AU 

Measurements of density fluctuations in the solar wind at 
~ 1 AU have been carried out by a variety of methods (e.g., 
Celnikier et al. 1983, 1987; Hnat et al. 2005; Kellogg & 
Horbury 2005); these allow us to calculate an upper limit on 
the heating by low-frequency Alfvenic turbulence analogous 
to that derived in the previous section. 

For concreteness, we use Celnikier et al. (1987)'s period II 
to find an upper limit to the KAW contribution to the density 
fluctuations, but similar results are obtained from other mea- 
surements. From their Figure 7, we infer that 8n^/«o ^ 10~ 2 
at k ± pt ~ 0.3, 12 where B ~ 1.5 x 10~ 4 G, no — 18 cm~ 3 , T p ~ 
1 .5 x 10 5 K, p p ~ 2 x 10 6 cm, v A ~ 77 km s" 1 , v wind ~460 km 
s _1 , and p ~ 0.43 in this epoch of data. Unl ike in the obser- 
vations of the solar corona described in 33.11 there is no clear 
evidence for an inner scale to the density fluctuations at 1 AU. 
Given the lack of a preferred scale, we evaluate the density 
fluctuations at k±pi ~ 0.3 as a compromise between the scales 
where the KAW density fluctuations peak (k±pi ~ 1) and the 
scales where incompressible MHD is a reasonable model for 
the cascade (k±pj <C 1); our conclusions are, however, insen- 
sitive to reasonable variations about this choice. 

A limit on the density fluctuations due to KAWs of 
Sn^/no < 10~ 2 at k±pi ~ 0.3 corresponds to a limit on the 
velocity fluctuations of Sv^/va < 10~ 2 at the same scale (i.e., 
§ v /t S 1 k m s _1 ) and thus to a limit on the heating rate due to 
low-frequency Alfvenic turbulence of 

e<5 x lO-^ergcnrV 1 . (13) 

1 ' For comparison, the heating rate at r = 5R . was 3 x 10~ 10 erg cm -3 s -1 
in the theoretical model of Markovskii & Hollweg (2002), and between 8 x 
10~ 10 ergcm~ 3 s~' and 5 x 10~ 9 ergcm~ 3 s~' in the "sweeping" models of 
Marsch & Tu (1997). 

12 This corresponds to a satellite-frame frequency ~ 1 Hz. 



At ~ 1 AU, the heating rate in the solar wind can be directly 
measured from the non-adiabatic temperature profile of the 
protons and electrons, using e ~ pv r Tds/dr. For example, us- 
ing the proton data from Voyager 2 (e.g., Matthaeus et al. 
1999), we infer that £ ~ 3 x 10~ 16 ergcm _3 s _1 is required 
to explain the non-adiabaticity of the solar wind at 1 AU. 13 
The close correspondence between this measured heating rate 
at 1 AU and the upper limit in equation dT3b implies that, if 
low frequency Alfvenic turbulence contributes to heating the 
solar wind at ~ 1 AU, direct density fluctuations due to com- 
pressive KAWs must contribute significantly to the measured 
density fluctuations at k±pj ~ 1. In this context, we note that 
Kellogg and Horbury (2005) have argued for an ion-acoustic 
or KAW origin for the small-scale density fluctuations in the 
solar wind at 1 AU using completely independent arguments. 
It is also important to note that the Celnikier et al. mea- 
surements are in the relatively slow solar wind (v ~ 450 km 
s ), in which the turbulence is observed to be fairly balanced 
(Grappin et al. 1990). Thus the cascade time is likely com- 
parable to the lower limit in equation dlPl l, in which case the 
cascade power can also be comparable to the upper limit in 
equation ( fT3l . 

The density fluctuation measurements of Celnikier et al. 
(1987) show a break from a Kolmogorov spectrum at low 
k to a flatter spectrum at high k. This is qualitatively con- 
sistent with the expectations from Figure [TJ The observed 
break happens at a rest frame frequency of ~ 0.1 Hz (Cel- 
nikier et al. 1987), which corresponds to k^ 1 ~ 7 x 10 cm 
3> Pi,di ~3x 10 6 cm using the Taylor hypothesis. This sug- 
gests that / <C 1, i.e., that the passive scalar contribution to 
the density fluctuations is small compared to the KAW contri- 
bution, so that the latter begins to dominate at k±pt <C 1 (see 
Fig. [TJ. A small value of / at 1 AU is not implausible given 
the observations described above, following equation (01. 

4. OBSERVATIONAL CONSTRAINTS ON THE "SWEEPING" MODEL 
OF CYCLOTRON HEATING 

An alternative coronal heating mechanism that has been 
considered in some detail is ion-cyclotron heating by high- 
frequency (kHz-range) waves. One of the proposed sources 
for these kHz-range waves is magnetic reconnection in the 
photosphere or chromosphere (Axford & McKenzie 1992). 
Because Bq and £2/ decrease with increasing r, waves with 
CO < £2, at the base of the corona eventually reach radii at 
which CO ~ £2,, at which point the waves undergo cyclotron 
damping and heat the ions. If a broad frequency spectrum of 
waves is launched from the Sun, the waves will result in radi- 
ally extended ion heating, with lower-frequency waves caus- 
ing heating farther from the Sun (Schwartz et al 198 1 ; Axford 
& McKenzie 1992; Marsch & Tu 1997; Tu & Marsch 1997; 
Ruzmaikin & Berger 1998; Czechowski et al 1998). Although 
cyclotron heating could in principle explain the observed tem- 
perature anisotropies of ions in the corona and fast solar wind, 
this "sweeping" or "direct-launching" model faces a signifi- 
cant difficulty. If the waves are oblique, with wave vectors 
that make a nonzero angle 9 with respect to Bo, then the 
waves induce density fluctuations at high frequencies. Holl- 

13 Cranmer & Van Ballegooijen 2005 and Cranmer et al (2009) find sim- 
ilar heating rates at 1 AU. Celnikier et al. (1987) do not report the electron 
temperature for their solar wind epochs, but the mean electron temperature 
for this solar wind speed is ~ 1.5 x 10 s K (Newbury et al. 1998), which is 
comparable to T p . We thus expect the electron heating rate to be at most com- 
parable to the proton heating rate; the proton heating rate is thus a reasonable 
proxy for the total heating rate, at the factor of 2 level. 
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weg (2000) investigated this effect, modeling the wave obliq- 
uity by setting 9 = 60°. Using the wave power spectra em- 
ployed in the "sweeping" models of Marsch & Tu (1997), he 
found that the waves would induce larger density fluctuations 
than are detected by radio observations, by a factor of > 10 2 
at k = 0.3 x 10~ 5 cm" 1 . Tu & Marsch (2000, 2001) have 
responded to this criticism by suggesting that transit-time 
damping could remove the obliquely propagating waves and 
thereby reduce the density fluctuations. We believe, however, 
that "phase mixing" by laminar (Heyvaerts & Priest 1983) and 
turbulent (Chandran 2008b) density structures will transfer 
Alfven-wave energy to larger k±, thereby significantly limit- 
ing the amount of Alfven/ion-cyclotron wave energy that can 
remain at small 9 as the waves propagate from the coronal 
base out to r = 5Rq. 

4.1. Why Radio Observations Rule Out the "Sweeping" 
Model but not Turbulent Heating 

In the "sweeping" or "direct-launching" model, the heating 
rate from high-frequency waves 2hf satisfies 



Air 



'pr 



(14) 



where Ehf is the energy density of waves with frequencies 
between 0.5£2; and £2, and t pr is the time for Alfven waves to 
propagate through the distance over which the magnetic field 
strength decreases by a factor of 2. Neglecting the solar-wind 
bulk-flow speed, and assuming Bq °c r~ 2 , the value of f pr is 
~ 0.4r/vA, which is ~ 800 s at r = 5R, . in our model coronal 
hole. On the other hand, in the turbulent-heating model, the 
heating rate for low-frequency turbulence 2if satisfies 



Gif ' 



Elf 

tc ' 



(15) 



where Ejf is the energy density of low-frequency Alfvenic 
fluctuations with 0.5A:,- < k± < ki, and ? c is the cascade time 
at k± = k(. Equations (0 and ( TTOb imply t c > 0.1 s, so 
that t c may be as much as four orders of magnitude smaller 
than f pr . For a fixed heating rate and ? c <C ? pr , the turbulent- 
heating model requires much less fluctuation energy at small 
scales than the sweeping model. In addition, because the 
compressibility of anisotropic KAWs is comparable to that 
of oblique ion-cyclotron waves, the turbulent-heating model 
also requires much lower density fluctuations than the sweep- 
ing model for a fixed heating rate when t c <C f pr - 

5. CONCLUSION 

In this paper, we have used density fluctuation measure- 
ments to test models in which the solar wind is heated by 
low-frequency Alfven-wave turbulence; the turbulent fluctu- 
ations are assumed to satisfy the inequality k± ^> \k\\\. The 
density fluctuation measurements - many of which are based 
on scintillation at radio wavelengths - place an upper limit on 
the amplitude of Alfven waves at small scales, where oblique- 
Alfven-wave compressions (kinetic Alfven waves = KAWs) 
induce density fluctuations. This upper limit in turn implies 
an upper limit on the turbulent heating rate. We calculate this 
upper limit for coronal holes at r — 5Rq and in the near-Earth 
solar wind at 1 AU. The upper limit on the turbulent heating 
rate we derive (eq. ffTTI ) can be realized only if two conditions 
are satisfied: (1) most of the measured small-scale density 
fluctuations are in fact from KAWs, and (2) the timescale for 
energy to cascade from one scale to another is ~ (£iSv;fc ± ) . 



The latter is only expected to be the case if there is compara- 
ble energy in Sunward and anti-Sunward waves in small-scale 
turbulent fluctuations. It is currently unclear whether this is in 
fact the case, both theoretically and empirically (especially 
at ~ 5Rq), which is one of the primary reasons that density 
fluctuations can only place an upper limit on the heating rate 
produced by low-frequency Alfvenic turbulence (in their re- 
lated analysis Coles & Harmon 2005 implicitly assumed that 
the turbulence is balanced). 

The upper limit on the turbulent heating rate at r = 5Rq 
in coronal holes exceeds the parameterized turbulent heating 
rates employed in models of the fast solar wind , by a factor of 
a few to ~ 100 depending on the models ( 33. It . At 1 AU, on 
the other hand, the upper limit on the turbulent heating rate 
due to low-frequency Alfvenic turbulence is within a factor 
of 2 of the measured heating rate - the latter being from the 
measured non-adiabatic temperature profiles ( 13.21 . We thus 
conclude that models in which the solar wind is accelerated by 
heating from low-frequency Alfvenic turbulence are consis- 
tent with the measured density fluctuations in the solar wind. 
By contrast, models in which the solar wind is accelerated by 
non-turbulent high frequency waves "sweeping" through the 
ion-cyclotron resonance are inconsistent with the measured 
density fluctuations at 5Rq unless the ion-cyclotron waves 
have k± = to very high precision (Hollweg 2000). The 
key difference between the turbulent model and the sweep- 
ing model is that in the sweeping model the timescale to dis- 
sipate the energy contained in small-scale (~ p,) fluctuations 
is set by the expansion speed of the solar wind, while in the 
turbulence mode l it i s set by the much faster cascade time 
at small scales ( 34. 11 1. Thus for a fixed heating rate the tur- 
bulent heating model has a much smaller energy density in 
small-scale fluctuations, which in turn produce much smaller 
density fluctuations. 

Future measurements and calculations could help to clar- 
ify the relative importance of passive fluctuations and KAW 
compressions in producing density fluctuations at different 
scales for different values of p. For example, the passive 
scalar contribution to the density fluctuations in the solar wind 
should be suppressed when p > 1, because the damping rate 
of the passive fluctuations (slow waves and entropy modes) 
is larger than the cascade rate (§0. Thus, measurements of 
the density fluctuations in solar wind epochs with large p 
are likely to be particularly instructive; more detailed theo- 
retical calculations of the suppression of the passive scalar 
contribution at p > 1 would aid in interpreting such measure- 
ments. In addition, future progress towards understanding the 
inertial-range power spectra of the density fluctuations, Sun- 
ward Alfven waves, and anti-Sunward Alfven waves in un- 
balanced (or cross-helical) turbulence could lead to signifi- 
cantly tighter constraints on the heating rate contributed by 
KAW turbulence, particularly near the Sun where the energy 
in anti-Sunward Alfven waves greatly exceeds the energy of 
Sunward Alfven waves. Finally, the fact that the upper limit 
on the KAW heating rate at 1 AU derived from density fluctu- 
ations is nearly equal to the measured heating rate in the solar 
wind ( 33.2b strongly motivates a more detailed analysis of the 
small-scale density fluctuations and their implications. 



This work was supported in part by the the Center for 
Integrated Computation and Analysis of Reconnection and 
Turbulence (CICART) under DOE under grant number DE- 
FG02-07-ER46372, by the NSF/DOE Partnership in Basic 



8 



Plasma Science and Engineering under grant number AST- 
0613622, and by NASA under grant numbers NNX07AP65G 
and NNH06ZDA001N-SHP06-0071 at the University of New 



Abraham-Shrauner, B. & Feldman, W. C. 1977, J. Geophys. Res., 82, 618 
Allen, L. A., Habbal, S. R., & Hu, Y. Q. 1998, J. Geophys. Res., 103, 6551 
Antonucci, E., Dodero, M. A., & Giordano, S. 2000, Solar Phys., 197, 1 15 
Axford, W. I. & McKenzie, J. F. 1992, in Solar Wind Seven, ed. Marsch, E. 
& Schwenn, R. R. (New York: Pergamon, p. 1) 
Barnes, A. 1966, Phys. Fluids, 9, 1483 

Bavassano, B., Pietropaolo, E., & Bruno, R. 2000, J. Geophys. Res., 105, 
15959 

Belcher, J. W., & Davis., L. J., & Smith, E. J. 1969, J. Geophys. Res., 74, 
2302 

Beresnyak, A., & Lazarian, A. 2008, ApJ, 682, 1070 
Boldyrev, S. 2006, Phys. Rev. Lett., 96 115002 

Celnikier, L. M., Harvey, C. C, Jegou, R., Kemp, M., & Moricet, P. 1983, 
A& A, 126, 293 

Celnikier, L. M., Muschietti, L., & Goldman, M.V. 1987, A& A, 181, 138 

Chandran, B. D. G. 2005, Phys. Rev. Lett., 95, 265004 

Chandran, B. D. G. 2008a, ApJ, 685, 646 

Chandran, B. D. G. 2008b, Phys. Rev. Lett., 101, 235004 

Chandran, B. D. G., Quataert, E., Howes, G. G., Hollweg, J. V., and Dorland, 

W. 2009, ApJ, in press. 

Coleman, P. J. 1968, ApJ, 153, 371 

Coles, W. A., & Harmon, J. K. 1989, ApJ, 337, 1023 

Coles, W. A., Grail, R. R., & Klinglesmith, M. T. 1995, J. Geophys. Res., 
100, 17069 

Cranmer, S. R. 2000, ApJ, 532, 1197 

Cranmer, S. R. & van Ballegooijen, A. A. 2003, ApJ, 594, 573 

Cranmer, S. R. & van Ballegooijen, A. A. 2005, ApJS, 156, 265 

Cranmer, S. R., van Ballegooijen, A. A., & Edgar, R. J. 2007, ApJS, 171, 520 

Cranmer, S. R., Matthaeus, W. H., Breech, B., & Kasper, J. C. 2009, ApJ, in 

press 

Czechowski, A., Ratkiewicz, R., McKenzie, J. E, & Axford, W. I. 1998, A& 
A, 335, 303 

Dmitruk, P., Milano, L. J., & Matthaeus, W. H. 2001, ApJ, 548, 482 
Dmitruk, P., Matthaeus, W. H., Milano, L. J., Oughton, S., Zank, G. P., & 
Mullan, D. J. 2002, ApJ, 575, 571 

Dmitruk, P., Matthaeus, W. H., & Seenu, N. 2004, ApJ, 617, 667 
Dobrowolny, M., Mangeney, A., Veltri, P. L. 1980, Phys. Rev. Lett., 76, 3534 
Feldman, W. C, Habbal, S. R., Hoogeveen, G., & Wang, Y.-M. 1997, J. 
Geophys. Res., 102, 26905 

Fisher, R. & Guhathakurta, M. 1995, ApJL, 447, 139 
Fisk, L. A. 2003, J. Geophys. Res., 108, 7 

Fisk, L. A., Gloeckler, G., Zurbuchen, T. H., Geiss, J., & Schwadron, N. A. 
2003, Solar Wind Ten, Am. Inst. Phys. Conf. Ser., 679, eds. Velli, M., Bruno, 
R., Malara, E, & Bucci, B., 287 

Gary, S. P. & Nishimura, K. 2004, J. Geophys. Res., 109, 2109 

Goldreich, P., & Sridhar, S. 1995, ApJ, 438, 763 

Grappin, R., Pouquet, A., & Leorat, J. 1983, A&A, 126, 51 

Grappin, R., Mangeney, A., & Marsch, E. 1990, J. Geophys. Res., 95 8197 

Gruzinov, A. 1998, ApJ, 501, 787 

Harmon, J., & Coles, W. 2005, J. Geophys. Res., 110, A03101 
Hartle, R. E., & Sturrock, P. A. 1968, ApJ, 151, 1155 

Hellinger, P., Travmcek, P., Kasper, J. C, & Lazarus, A. J. 2006, 

Geophys. Res. Lett., 33, 9101 

Heyvaerts, J. & Priest, E. R. 1983, A& A, 1 17, 220 

Hollweg, J. V. 1999, J. Geophys. Res., 104, 14811 

Hollweg, J. V. 2000, J. Geophys. Res., 105, 7573 

Hollweg, J. V., & Isenberg, P. A. 2002, J. Geophys. Res., 107, 1 

Hollweg, J. V. & Turner, J. M. 1978, J. Geophys. Res., 83, 97 

Holzer, T. 1977, J. Geophys. Res., 82, 23 

Hossain, M., Gray, P. C, Pontius, D. H., Matthaeus, W. H., & Oughton, S. 
1995, Phys. Fluids, 7, 2886 

Howes, G. G., Cowley, S. C, Dorland, W., Hammett, G. W., Quataert, E., & 
Schekochihin, A. A. 2008a, J. Geophys. Res., 113, A05103 
Howes, G. G., Dorland, W., Cowley, S. C, Hammett, G. W., Quataert, E., 
Schekochihin, A. A., & Tatsuno, T. 2008b, Phys. Rev. Let., 100, 65004 



Hampshire. E. Quataert was supported in part by NSF-DOE 
Grant PHY-081281 1 and by NSF Grant ATM-0752503. 



Iroshnikov, P. 1963, Astron. Zh. 40, 742 

Isenberg, P. A. & Hollweg, J. V. 1983, J. Geophys. Res., 88, 3923 
Johnson, J. R. & Cheng, C. Z. 2001, Geophys. Res. Lett., 28, 4421 
Klinglesmith, M. 1997, "The Polar Solar Wind from 2.5 to 40 Solar Radii: 
Results of Intensity Scintillation Measurements," PhD Thesis, Univ.of Calif., 
San Diego 

Kohl, J. L. et al 1998, ApJL, 501, 127 
Kraichnan, R. H. 1965, Phys. Fluids, 8, 1385 

Leamon, R. J., Smith, C. W., Ness, N. E, & Wong, H. K. 1999, J. Geophys. 
Res., 104, 22331 

Lithwick, Y, & Goldreich, P. 2001, ApJ, 562, 279 
Lithwick, Y, & Goldreich, P. 2003, ApJ, 582, 1220 
Lithwick, Y, Goldreich, P., & Sridhar, S. 2007, ApJ, 655, 269 
Little, L. T., & Ekers, R. D. 1971, A&A, 10, 306 
Markovskii, S. A. 2004, ApJ, 609, 1112 

Markovskii, S. A., & Hollweg, J. V. 2002a, Geophys. Res. Lett., 29, 24 
Markovskii, S. A., & Hollweg, J. V. 2002b, J. Geophys. Res., 107, 21 
Markovskii, S. A., Vasquez, B. J., Smith, C. W., & Hollweg, J. V. 2006, ApJ, 
639, 1177 

Marsch, E. 1991, in Physics of the Inner Heliosphere, ed. Schwen, R. & 
Marsch, E. (New York: Springer- Verlag: p. 159) 
Marsch, E., & Tu, C.-Y. 1997, A&A, 319, L17 

Marsch, E., Ao, X.-Z., & Tu, C.-Y. 2004, J. Geophys. Res., 109, 4102 
Maron, J. & Goldreich, P. 2001, ApJ, 554, 1175 

Matthaeus, W. H., Zank, G. P., Leamon, R. J., Smith, C. W., Mullan D. J., & 
Oughton, S. 1999, Sp. Sci. Rev., 87, 269 

Matthaeus, W. H., Zank, G. P., Smith, C. W., & Oughton, S. 1999, Physical 
Review Letters, 82, 3444 

Matthaeus, W. H., Mullan, D. J., Dmitruk, P., Milano, L., & Oughton, S. 2003, 
Nonl. Proc. Geophys., 10, 93 

McKenzie, J. F, Ip, W. H., & Axford, W. I. 1979, Astrophys. Sp. Sci., 64, 
183 

Newbury, J. A., Russell, C. T, Phillips, J. L., & Gary, S. P. 1998, 

J. Geophys. Res., 103, 9553 

Parker, E. 1965, Space Sci. Rev., 4, 666 

Perez, J., & Boldyrev, S. 2009, Phys. Rev. Lett., 102, 025003 

Podesta, J., & Bhattacharjee, A., Phys. Rev. Lett., submitted 

Quataert, E. 1998, ApJ, 500, 978 

Quataert, E., & Gruzinov, A. 1999, ApJ, 520, 248 

Roberts, D. A., Goldstein, M. L., Klein, L. W., & Matthaeus, W. H. 1987, J. 
Geophys. Res., 92, 12023 

Ruzmaikin, A. & Berger, M. A. 1998, A&A, 337, L9 

Schekochihin, A. A., Cowley, S. C, Dorland, W., Hammett, G. W., Howes, 
G. G., Quataert, E., & Tatsuno, T. 2009, ApJS, 182, 310 
Sahraoui, E, Goldstein, M. L., Robert, P., & Khotyaintsev, Y. V. 2009, 
Physical Review Letters, 102, 231102 

Schwartz, S. J., W. C. Feldman, & Gary, S. P. 1981, J. Geophys. Res., 86, 541 
Scudder, J. D. 1992a, ApJ, 398, 299 
Scudder, J. D. 1992b, ApJ, 398, 319 

Shebalin, J. , Matthaeus, W., & Montgomery, D. 1983, J. Plasma Phys., 29, 
525 

Tu, C. Y, & Marsch, E., 1990, J. Geophys. Res., 95, 4337 

Tu, C. Y, & Marsch, E., 1995, Sp. Sci. Rev., 73, 1 

Tu, C.-Y, & Marsch, E. 1997, Solar Phys., 171, 363 

Tu, C.-Y, & Marsch, E. 2001a, J. Geophys. Res., 106, 8233 

Tu, C.-Y, & Marsch, E. 2001b, A&A, 368, 1071 

Tu, C.-Y, Pu, Z.-Y, & Wei, F.-S. 1984, J. Geophys. Res., 89, 9695 

Verdini, A., & Velli, M. 2007, ApJ, 662, 669 

Velli, M. 1993, A&A, 270, 304 

Voitenko, Y. & Goossens, M. 2004, ApJL, 2004, 605, 149 
Zhou, Y, & Matthaeus, W. H. 1990, J. Geophys. Res., 95, 14881 



